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ABSTRACT 

3 
The  He      (n,  p)T  reaction  can  be  utilized  in  a  proportional  counter 

for  neutron-energy  determination  in  the  region   100  kev  to    1   Mev.      This 

has  already  been  demonstrated  and  the  reaction  cross   section  in  this 

energy  region  has  been  measured  with  a  counter.      For  practical 

application  of  the  counter  as  a  spectrometer,    such  as  in  health  physics 

work,    it  is  desirable  to  increase  the  counter  efficiency.      A  proportional 

counter  utilizing  an  anticoincidence  ring  to  reduce  wall  effect  has  been 

constructed  for  this  purpose.      Some  preliminary  tests  of  the  proposed 

system  are  described. 
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CHAPTER  I 

INTRODUCTION 

The  detection  of  neutrons  and  especially  the  measurement  of 
their  energy  have  from  the  first  attempts  presented  a  more  difficult 
problem  than  for  the  other  common  nuclear  particles.      This  is  be- 
cause the  usual  methods  depended  on  effects  resulting  from  the  charges 
of  the  particles.      The  experiments  by  Dee  in  which  he  investigated  the 
ionization  produced  in  air  in  a  cloud  chamber  irradiated  by  fast 
neutrons'-        were  reported  at  the  same  time  as  Chadwick's  announce- 
ment of  the  discovery  of  the  neutron.      Dee  concluded  that,    if  the 
neutron  interacts  with  atomic  electrons  at  all,    the  process  produces 
not  more  than  one  ion  pair  per  three  meters  of  the  neutron's  path. 

If  we  accept  an  arbitrary  description  of  neutron  energy,    E    ,    as 
intermediate  neutrons:    1  kev  to  500  kev,    fast  neutrons:   0.5  Mev  to 
10  Mev,    then  we  may  say  in  general  that  the  predominant  reaction  of 
intermediate  neutrons  with  nuclei  is  elastic   scattering,    and  that  in  the 
fast  range  many  other  reactions  appear,    the  most  important  of  which 
is  inelastic  scattering1       .      Hydrogen-   or  methane-filled  ionization 
chambers  or  proportional  counters  became  an  important  method  of 
detecting  fast  neutrons  because  the  neutron  can  impart  practically  ail 
of  its  energy  to  a  proton  in  a  head-on  collision,    and  the  recoil  proton 
is  then  the  particle  causing  the  action  of  the  counter.     However,    a 
determination  of  the  neutron  energy  spectrum  by  the  recoil  method 
requires  a  thin  radiator  and  a  double  collimation- -first  of  the  neutron 
beam  and  then  of  the  recoil  protons- -  resulting  often  in  too  low  an 
efficiency.      Theoretically,    the   second  collimation,    may  not  be 
necessary  because  the  neutron  spectrum  can  be  obtained  from  the 
recoil  spectrum  by  differentiation,    but  this  leads  to  large  errors  since 
we  are  effectively  taking  the  differences  of  numbers  of  about,  the  same 
size.     Nuclear  emulsions  and  organic  scintillation  counters  have  be- 
come very  important  detectors  in  neutron  spectroscopy,    especially  for 
neutron  energies  above   1  Mev, 
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For  energies  below  1  Mev  it  became  apparent  that  exothermic 
nuclear  reactions  might  be  useful  for  neutron  spectroscopy,      New 
methods  were  especially  needed  in  health  physics  applications,    where 
any  schemes  requiring  collimation  are  automatically  ruled  out.      Fast 
neutrons  are  of  major  concern  to  health  physicists  because  of  the 

difficulty  of  shielding  against  them,    and  their  large  relative  biological 

[3] 
effectiveness        .     In  1946  Feld  proposed  that  a  proportional  counter 

utilizing  the  (n,  p)  reaction  might  be  useful  for  neutron  spectroscopy 

in  the   10-  to   1000-kev  range     with  the  advantage  over  recoils  of 

Ul 
detecting  lower  energies  and  making  collimation  unnecessary        .      He 

considered  the  most  promising  reactions  at  that  time: 

14  14 

CI        (n9  p)S 

However,    as  he  expected,    these  do  not  very  well  fulfill  the  first  of 
the  general  requirements,    which  might  be  listed  as  follows: 

(1)  The  reaction  cross  section  should  be  fairly  large  and  its 
variation  with  neutron  energy  smooth  and  free  from  resonances  in 
the  energy  region  of  interest  (and  accurately  known). 

(2)  The  Q  value  of  the  reaction  should  be   small,    because   (E      +  Q) 
is  the  total  energy  measured  by  a  counter,    and  for  accuracy  E 
cannot  be  too  small  a  fraction  of  the  total. 

(3)  There   should  be  no  low-lying  excited  states  of  the  residual 
nucleus  to  cause  an  ambiguity  in  the  energy  of  the  reaction  particles, 

We  can  look  for  these  requirements  to  be  met  over  appreciable 
ranges  only  in  light  nuclei,    in  which  the  energy  levels  are  more 
widely  spaced,    avoiding  resonances  and  excited  states.      Also,    for 
small  Q  values  appreciable  penetration  of  the  Coulomb  barrier  for 
(n,  p)  or   (n,  a)  charged-particle  reactions  would  be  limited  to  light 
nuclei. 

Two   (n,  a)  possibilities  are 

Ei6   (n,  a)T, 

B        (n,  a)Li    . 
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The  first  has  a  large  Q  (4.78  Mev)  and  a  resonance  at  E     -  0.27  Mev. 
The  second,    widely  used  for  thermal  neutron  detection,    has  an  ex- 
cited state  and  fairly  large  Q  values   '2.79  and  2.31  Mev). 

With  supplies  of  helium- 3  becoming  more  available  and  a  likely 
3 
possibility  being  He     n,  p)T,    helium- 3-filled  proportional  counters 

came  under  consideration.      In  1950  Coon  showed  that  the  cross 

[Si 
section  in  the   1-Mev  region  was  about   1  barn        ,    comparable  to  the 

scattering  cross   section  of  hydrogen.      The  thermal  cross   section  for 

[61 
this  reaction  is  large   (5400  ±  300  barns)        .      These  types  of  reaction 

cross  sections  obey  the   l/v  law  for  a  wide  range,    but  intermediate 

neutrons  may  depart  from  the   l/v  law. 

Some  work  on  helium- 3-filled  proportional  counters  has  been 

[7l 
done  by  Batchelor  and  others  at  Harwell,    England1       .      They  used  a 

small  amount  of  helium- 3  admixed  with  xenon  or  krypton  as  the 

counting  gas  and  stopping  gas  to  reduce  wall  effect.      Also,    a  small 

amount  of  carbon  dioxide  was  added  to  stabilize  the  multiplication 

process.      In  their  final  filling  the  pressure  of  krypton  used  was   164 

cm  of  mercury.      Since  they  found  that,    for  E     =   1  Mev,    a  pressure 

of  8  atmos  of  krypton  is  required  to  reduce  the  wall  effect  to   10%, 

they  set  up  a  computer  program  to  correct  for  wall  effect.      Their 

work  remeasured  the  neutron  cross  section  in  the  range   120  kev  to 

1  Mev  more  accurately  than  previously,    showing  a  much  flatter 

energy  dependence  than  the    l/v  law. 

They  also  pointed  out  a  basic  source  of  ambiguity  in  energy 

determination  other  than  counter  defects       As  we   show  in  Chaper  VIII, 

the  maximum  recoil  energy  of  a  He     nucleus  is  ER    ^  ~  3/4  E    .      The 

true  energy  counts  run  from  Q  for  E     -  0  to  E     +  Q.      Thus,    for 

3/4  E     =  Q,    the  recoil  spectrum  begins  to  overlap  the  reaction 

spectrum.      On  the  basis  of  Batchelor's  value  of  Q  as  770  kev,    this 

means  that  ambiguities  exist  if  neutrons  of  greater  energy  than  1.03 

Mev  are  present.      The  scattering  cross  section  up  to  about  20  Mev  is 

comparable  in  size  to  the  reaction  cross  section  of  interest,    but  of 

course  only  a  fraction  of  those  neutrons  scattered  deliver  a  large 
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3 
portion  of  their  energy  to  the  He      recoils.      In  any  case,    the  presence 

of  high-energy  neutrons  reduces  the  accuracy. 

Gaseous  counters  have  notoriously  low  efficiency,    and  as 
helium-3  has  become  more  plentiful  and  less  expensive,    it  seems 
worthwhile  to  investigate  the  feasibility  of  increasing  the  counting 
rate  by  increasing  the  quantity  of  helium-3  in  the  counter.      As  we  go 
to  higher  pressures,    special  problems  arise,    such  as  the  demands 
for  gas  purity,    higher  voltages,    and  better  insulation1-       .      Also,    the 
container  must  be  more  massive,    and  proper  consideration  should 
be  given  to  the  effects  of  the  masses  upon  the  spectral  distribution 
of  the  neutrons,    primarily  caused  by  inelastic  collisions   (trans- 
parency correction).       For  these  reasons  it  would  be  desirable  to 
use  helium~3  alone  as  the  primary  filling,    utilizing  it  for  the  count- 
ing and  stopping  action.      However,    we   soon  find  it  very  difficult  to 
decrease  wall  effect  to  a  reasonable  value  by  helium-3  pressure 
alone  because  stopping  powers  decrease  with  lower  atomic  numbers. 

For  instance,    the  following  ranges  of  a  5.3-Mev  alpha  particle  at  one 

[9I 
atmosphere  have  recently  been  measured        : 

helium- -Zl  cm, 

krypton- -3.0  cm, 

xenon-  -  2.  2  cm. 
For  this  reason  it  is  desirable  at  the   same  time  to  determine  the 
practicality  of  reducing  the  wall  effect  by  surrounding  the   sensitive 
volume  of  the  counter  with  an  anticoincidence  ring. 

Since  helium- 3  is  quite  valuable,    a  fairly  elaborate  filling 
system  must  be  devised  to  insure  no  loss  of  gas  in  filling  and  sub- 
sequent recovery.      The  purity  requirements  are  not  so  great  as  in 
noble  gas   scintillation  counting,    in  which  a  small  amount  of  contam- 
inant such  as  oxygen  or  other  polyatomic  gases  may  completely 
quench  the  useful  scintillation  property,    but  some  purification  arrange 
ments  will  probably  be  needed  in  the  filling  system  to  achieve  good 
energy  resolution  and  to  remove  impurities  which  inevitably  get  in 
to  the  gas  during  filling  and  recovery.      It  is  important  that  the 


loss  of  helium-3  be  kept  small  during  the  purification  process.      Purifi- 
cation procedures  for  noble  gases  using  metallic  calcium  have  been 
applied  successfully  for  some  time  in  removing  all  electroi-negative 
impurities,    nitrogen,    hydrogen,    and  to  a  lesser  degree  carbon  dioxide 
Also,    the  method  of  physical  adsorption  of  impurities  in  an  activated 
charcoal  trap  cooled  to  liquid  nitrogen  temperatures  has  been  used 

successfully  in  a  helium-3  system,    and  has  the  advantage  of  greater 

[111 
ease  and  convenience1         .      However,    it  appears  there  maybe  a  special 

problem  due  to  the  radioactive  background  from  tritium  contamination 

of  helium-3,    and  this  problem  is  increased  by  using  higher  helium-3 

pressure 

It  was  felt  that  useful  information  about  these  two  techniques,    (a) 

using  helium-3  as  the  detecting  medium  and  primary  counting  gas  and 

b     reducing  wall  effect  with  an  anticoincidence  ring,    could  be  obtained 

by  testing  a  counter  with  a  filling  of  ordinary  helium.      The  natural 

abundance  of  helium-3  is  about  0.00013,    so  that  it  is  not  ordinarily 

noticeable.      Owing  to  the  large  cross  section  at  thermal  energies,    it 

can  be  detected  in  tank  helium  by  use  of  a  nuclear  reactor  as  a  source 

of  neutrons,    and  this  has  been  suggested  as  a  way  of  surveying  world 

[12] 
helium  supplies  for  heium-3  content  .      Because  it  is  a  double- 

4 
magic-number  nucleus,    He     is  very  stable,    and  the  detection  of 

neutrons  in  a  helium-filled  counter  is  due  to  the  recoil  alpha  particles. 

Helium-3  and  helium-4  have  the   same  electronic  configuration,    and 

we  can  obtain  information  about  the  counting  characteristics  of  one 

by  experimenting  with  the  other. 
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CHAPTER    II 

COUNTER  SYSTEM 

A  method  of  reducing  wall  effect  in  proportional  counters  is  to 
detect  the  particles  that  escape  from  the    sensitive  volume  and  exclude 
them  from  the  measurements.      To  accomplish  this  it  is  desirable  that 
the  walls  of  the  main  counter  be  transparent  and  completely  surrounded 
by  transparent-walled  counters  in  anticoincidence       The  feasibility  of 
Geiger-Mueller  counters  with  transparent  walls  --so  that  wall  thick- 
ness would  not  limit  the  radiation  to  that  above  a  minimum--   was 

investigated  in  1944,    and  it  was  found  that  symmetrically  placed   wires 

[131 
could  approximate  a  cylinder1  .      Single  counters  of  triangular  type 

were  constructed  with  good  results.      However,    when  two  such  counters 

were  placed  close  together,    there  was  a  deleterious  effect  due  to 

distortion  of  the  field  at  the  central  wire.      Better  results  were  obtained 

with  a  regular  hexagonal  shape. 

COUNTER  CONSTRUCTION 

A  counter  has  been  constructed  in  which  a  layer  of  gas  surround- 
ing the  main  counter  is  in  effect  turned  into  a  separate  detecting  device 
by  a  circular  array  of  wires.      The  geometry  is  shown  in  Figs.    1  and  2. 
The  inner  circle  of  wires,    grounded  to  the  cylinder,    serves    as  the  cathode 
for  the  main  counter  and  as  the  inner  cathode  wall  for  the  ring  counter. 
In  the  outer  circle  of  wires  alternate  wires   serve  as  cathode  and 
counting  wires  for  the  ring  counter.      Since  there  are  36  wires  in  each 
circle,    we  effectively  have  the  main  counter  surrounded  by   18  separate 
counter  s. 

An  arrangement  similar  to  this,    referred  to  as  a    "wall-less 

counter",    has  been  applied  to  low-energy  beta  spectroscopy  and  proved 

[141 
to  be  a  practical  instrument1         .      As  pointed  out,    this  prevents 

escaping  high-energy  particles  from  having  an  effect  on  the  low-energy 

end  of  the  spectrum,    but  the  effective   sensitive  volume  now  depends 

on  energy.      However,    inaccuracies  because  of  this  circumstance  will 
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Fig.     1.      Arrangement  of  counter, 
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Fig.    2.      Cross  section  of  counter  wiring. 
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be  much  less  than  from  uncorrected  wall  effect.      Another  application 
of  this  type  of  counter  geometry  is  the  so-called    "Giles  telescope" 

r  15] 

neutron  spectrometer  In  this  arrangement  the  counter  is  long 

and  narrow,    and  neutrons  are  collimated  along  the  long  axis.      The 
requirement  of  a  thin  radiator  is  obviated  by  using  gas  recoils  and  the 
anticoincidence  method  to  prevent  registration  of  pulses  that  recoil 
oblique  to  the  chamber  axis.      Monoenergetic  neutrons  then  have  a 
pulse-size  distribution  with  a  well-defined  peak,      A  recent  applica- 
tion of  this   scheme,    utilizing  hydrogen  and  propane  fillings  at 

various  pressures,    has  resulted  in  energy  spreads  of  10%  for  neutron 

T 16 1 
energies  from  a  few  hundred  kev  to   10  Mev 

Some  of  the  constructional  details  of  our  counter  follow.      The 
counter  cylinder  has   l/ 16-inch-thick  brass  walls  and  is  4  inches  in 
outside  diameter       The  total  gas  volume  less  the  inlet  piping  is   1.93 
liters.      The  vessel  has  been  pressure-  and  leak-tested  with  helium  at 
20  atmospheres,    making  operating  pressures  up  to   10  atmospheres 
possible  within  safety  requirements.      All  wires  are   stainless   steel  of 
3-mil  diameter.      The  central  wire  is  supported  by  kovar  insulators 
and  is  electrically  available  only  at  the  base.       Its  length  between 
thickened  supports  is  8.5  inches.      If  we  assume  that,    owing  to  a  fall- 
off  in  field  strength  near  the  ends,    only  90%  of  this  length  defines  the 
sensitive  volume,    the  sensitive  volume  of  the  main  counter  is  888  cc. 
The  ring-counter  anode  wires  and  all  cathode  wires  are  each  strung 
from  a  continuous  piece  of  stainless  steel  wire  and  threaded  through 
holes  in    ^-inch-thick  lava  insulators  at  top  and  bottom,    the  anode 
wire   system  being  electrically  available  through  another  kovar  at  the 
base,      A  major  drawback  to  this  method  of  stringing  is  the  difficulty 
of  repair  in  case  of  wire  breakage.      The  volume  defined  by  the  ring 
counter  is  600  cc. 

The  anticoincidence  ring  also  cancels  out  particles  that  enter 
the  counter  from  the  outside  or  are  emitted  owing  to  contamination  of 
the  brass  walls,      Korff  says  a  well-cleaned  surface   should  emit  0.36 

2  r  i7i 

alpha  particle     per  cm     per  hour  within  a  factor  of  ten1         ,    and 
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Sharpe  gives  0.09  alpha  particle  per  cm     per  hour  for  commercial 

f  18] 
brass  .      In  normal  operation  the  counter  will  probably  not  be 

sensitive  to  single  secondary  electrons  knocked  from  the  walls  by 

gamma  rays,    but  unwanted  counts  due  to  pile-up  will  be  more  likely 

in  the  ring  counter,    where  the  density  of  secondary  electrons     will  be 

greater,      A  proportional  counter  with  walls  of  copper  tubing   l/32 

inch  thick  has  been  used  for  counting  gamma  rays  of  energy  up  to   1.6 

[191 
Mev  ,      The  gamma-ray  attenuation      "tenfolding  length"  is 

approximately  60  g/cm     for  all  materials  in  the  energy  region  where 
Compton  scattering  is  the  most  important  process^         .      This  extends 
from  less  than  one  to  several  Mev  for  all  elements.      Multiplying  by 
density,    we  find  our  brass  wall  is   1.34  g/cm     thick.      The  tenfolding 
length  converted  to  a  l/e  length  is  26.1  g/cm    .      We  have 
exp  (-1.34/26,1)  =  95%,    therefore  about  5%  of  all  gamma  rays  imping- 
ing on  the  counter  wall  will  be  converted  to  secondary  electrons.      Of 
course,    many  of  these  electrons  will  be  absorbed  in  the  wall  before 
reaching  the  counter  gas,    but  gamma-ray  counting  efficiency  would 
be  reasonably  high  if  the  counter  were  operated  at  high  gas 
multiplication. 

The  voltage  at  which  a  counter  under  given  conditions  is  operated 

determines  the  gas  multiplication.      For  a  given  counter  Rossi  and 

[2ll 
Staub  have   shown  the  functional  relationship  of  multiplication  to  be 


,      Pa      , 


where 

a  =  wire  radius, 

b  =  cathode  radius, 

P  =gas  pressure, 

V  =»voltage  across  counter. 
Since  multiplication  usually  takes  place  in  a  cylindrical  region  of 
small  diameter  near  the  wire,    we  might  assume  our  two  counter 
systems  identical  for  purposes  of  figuring  multiplication,    except  for 
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the  difference  in  In  b/a.      The  multiwire  potential  problem  would  be 
very  difficult  to  compute,    but  if  we  take  an  average  b  value  for  one 
of  the  ring- counter  wires  (averaged  over  measurements  to  the  inner 
circle  of  wires,    to  an  adjacent  wire  in  same  ring     or  to  the  wall), 
for  the  ring  counter  we  obtain  b/a  =  225  and  for  main  counter  b/a  =   1000 
The  ratio  of  the  logarithms  is   1.27,    therefore  for  the  same  value 
of  gas  multiplication,    this  indicates  that  the  main  counter  should   have 
a  higher  voltage  by  a  factor  of  1.27.      Of  course,    the  determining 
factor  on  multiplication  of  the  ring  counter  will  be  to  make  it  as   sensi- 
tive as  possible  in  order  to  detect  all  particles  that  leave  the  main 
counter,    but  not  so  sensitive  that  a  large  number  of  background  counts 
keep  the  desired  counts  from  the  main  counter  gated   out. 

ELECTRONIC  CIRCUITRY 

A  block  diagram  of  the  electronic  circuitry  used  with  the 
counter  is  given  in  Fig.    3.      The  blocks  below  the  dotted  line  complete 
the  neutron- spectrometer  system,    but  investigation  of  the  operation 
of  the   system  was  made  without  this  by  scaling  both  coincidences  and 
anticoincidences.      A  list  of  units  actually  used  is  given  in  Table  I. 
One  channel  of  the  dual- channel  variable  gate  and  delay  unit  can  be 
used. for  the  variable-gate  block,    and  the  use  of  this  unit  for  the 
pulse-shaping  function  could  be  replaced  by  simpler  electronic  cir- 
cuitry.     This  unit  is  very  useful,    however,    for  experimenting  with 
different  values  of  discriminator  setting  and  of  delay  and  length  of 
output  pulses.      If  two  pulses  arrive  at  the  pulse- shaping  block  at 
about  the  same  time,    we  want  to  insure  that  there  is  no  anticoincidence 
pulse  to  gate  the  pulse-height  analyzer  on.      To  do  this  we  delay 
channel  No.     1   slightly  so  that  the  channel  No.    2  output  pulse  will 
definitely  be  at  the  coincidence  unit  before  the  channel  No.    1  pulse 
arrives.-    Also,    channel  No.    2  pulse  is  made  longer   so  that  the 
coincidence  unit  will  see  it  during  the  entire  time  that  it  sees  the 
channel  No.    1  pulse.      The  count  on  Scaler   1+2  indicates  the  number 
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Fig.    3.     Block  diagram  of  counter  electronic  circuitry. 
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Table  I.      Units  used  in  counter  electronic  circuitry. 


Function  Unit 


Regulated  High  Voltage  Supply  Northern  Scientific  Co.  ,    ser.     #  1 

Meter  and  Divider  Panel  UCRL  3  kv  Divider  Panel 

Linear  Amplifier,    ch.    #   1  UCRL  Mod  5,    Dwg.    1X5334 

Preamplifier,    ch.    #  1  UCRL  Dwg.    IX4353D 

Linear  Amplifier,    ch.    #  2  UCRL  Mod.    3,    Dwg.    2T4404 

Preamplifier,    ch.    #  2  UCRL  Dwg.    3T2393 

Pulse  Shaping  UCRL  Variable   Gate  and  Delay  Unit 

Dwg.    2T8084F  and  2T8154 

Coincidence   Unit  UCRL  8  Channel  Quadruple  Mixer 

Dwg.    3T5894-1 

Scaler   1  +  2  UCRL  1024  Scaler,    Mod.    2 

Dwg.    3T8934 

Scaler   1,    2  UCRL  2  Channel   1024  Scaler,    Mod.  1 

Dwg.    4T6895 
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of  pulses  from  channel  No,     1  not  registered  by  the  pulse-height 
analyzer.      The  versatile  variable  gate  and  delay  unit,    used  with  one 
input  pulse    and  in  conjunction  with  the  coincidence  unit.,    can  also  be 
used  as  a  single-channel  pulse-height  analyzer  if  desired. 

IONIC  MOTION 

With  this  size  of  counter  one  may  have  to  worry  about  the  time 
lag  between  formation  of  the  pulse  by  the  main  counter  and  by  the  ring 
counter  when  ionization  occurs  close  to  the  edge  of  the  main  counter. 
We  can  compute  the  travel  time  for  ions  if  we  know  the  mobility    and 
how  drift  velocity  depends  on  mobility  and  the  field.      Sharpe  has 

developed  convenient  equations  for  computing  ionic  motion  in 

f  18] 

chambers  .      In  general,    in  cylindrical  or   spherical  fields.,    the  field 


is 


where 


Vr 
e    - 


In  y 


y  =  b/a 
r  ~  radius, 

V"  voltage  across  counter, 
and  for  a  cylinder  n  =   1. 
An  ion  crosses  the  system  from  one  electrode  to  the    other  in  a  time 
given  by 

T  =  [anp+1  (ynp+1-l)]   /  [(np+1)  K  (V/ln  y)P]    , 


where  p  =    \  for  electrons  and  1  for  ions  in  the  equation  relating  drift 


velocity   w  to  mobility  K, 


w  -  KeP 


We  are  concerned  with  the  time  for  an  electron  to  travel  across  the 
counter  because  most  of  the  pulse  is  formed  near  the  wire   owing  to  the 
concentration  of  the  field  at  the  wire  whether  there  is  multiplication  or 
operation  as  an  ion  chamber.      If  we  take  a  case  of  V  -   1000  volts, 
P  =  2.5  atmospheres  of  helium,    and  use  Sharpe's  value  of 


-14- 


3/2 


K    =  3.5  •    10       cm   '     /  N/volt  -   secrVatmos  ,    we  find  for  electrons 
T  =  7.55    |j.sec.      However,    we  expect  that  the  mobility  of  a  helium- 
carbon  dioxide  mixture    may  easily  be  greater  by  a  factor  of  ten, 
based  on  the  way  CO,  increases  the  drift  velocity  of  electrons  in  other 
noble  gases.      This  is  explained  for  argon  on  the  basis  of  the  Ramsauer 
effect  and  the  high  first  excitation  potential  for  argon,     11.5  volts,    in 
comparison  with  the  low  excitation  levels  in  a  CO?  molecule  .      The 

first  excitation  potential  of  helium  is  even  higher,     19.77  volts, 
Extensive  work  has  been  done  on  collection  times  in  argon,    krypton, 
and  xenon,    and  it  is  reported  that  once  a  fraction  of  one  percent  of 

CO,  is  added,    the  effect  of  original  polyatomic  impurities  in  a  noble 

[22l 
gas  may  be  assumed  insignificant  .      The  transit  time  computed 

above  is  reduced  by  the  same  factor  as  that  by  which  the  mobility  is 

increased,    and  should  not  be  too  long  for  the  coincidence  requirements. 

The  transit  time  for  a  He     ion  in  pure  helium  to  cross  the  system 

from  wire  to  wall,    computed  for  the  same  conditions  as  for  the  electron 

and  on  the  basis  of  a  mobility  of  10.8  cm.     /volt- sec-atm.    is    1850  usee. 

This  time  is  not  too  important,    since  the  pulse  can  be  shortened  by 

shortening  the  amplifier  time  constant. 
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CHAPTER  III 
ENERGY  RESOLUTION 

A  discussion  of  energy  resolution  depends  on  the  complete  theory 

of  proportional- counter  action.      This  has  been  well  covered  in  several 

.'',  j       *••  i      [10,17,18,21.23,24,25,26]       •"  t, 

books  and  articles1  J  .      The  attempt  here  is 

to  point  out  important  factors  concerning  the  energy  resolution  of  this 

particular  counter,    using  developments  from  any  of  the  basic  sources 

but  modifying  the  terminology  as  necessary  in  order  to  be  consistent. 

INHERENT  LIMITATIONS 

Inherent  fluctuations  in  a  proportional  counter  are  from  (a)  the 
number  of  ion  pairs  released  by  monokinetic  radiation,    and  (b)  the 
size  of  the  avalanche  each  ionization  electron  produces.      From 
theoretical  work  based  on  the  assumption  that  the  first  was  a  Poisson 
distribution  (it  was  found  experimentally  to  be   smaller  than  this)  and 
expanded  upon  experimentally,    it  was  shown  that  relative  variance  of 
total  size  of  avalanche   is  approximately  given  by 

Vp/(mAm0)2     *     l/m0  (  HI- 1 ) 

where 

V     =  variance   (mean  square  deviation)  of  size  of  avalanche, 

m„  ~  mean  of  number  of  ions  initially  liberated, 

m  .  -  mean  of  number  of  ions  produced  in  avalanche  by  one 
electron. 

We  note  that  this  inherent  spread  is  no  longer  a  significant 

factor  at  the  energies  at  which  we  will  be  working  because  in  the   1-Mev 

energy  range  and  at  an  assumed  mean  energy  to  produce  an  ion  pair 

4 
of  30  ev,    m^  is  greater  than  10    .      These  theories  were  useful,    and 

experiments  were  done  down  in  the  kev  region  where  m0  was  a  smaller 
number.      For  our  magnitude  of  rn^,    the  relative  variance  would  be 
around  0.01%.      Since  experience  has  shown  that  the  relative  spread  of 
pulse  height  (relative  standard  deviation  of  pulse  height  from  a  mono- 
kinetic  radiation)  is  at  least   1%  or  2%  (owing  to  mechanical  defects  of 
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a  counter  or  defects  in  gas  falling)  at  these  energies,    the  inherent 
spread  is  not  a  significant  factor.      The  mean  energy  per  ion  pair, 
which  determines  m~,    is  expected  to  be  practically  constant  over  the 
(E      +  O)  range.      This  may  not  be  true  in  very  pure  helium,    in  which 
the  mean  energy  has  been  measured  as  41.3  ev,    in  contrast  to  29.7  ev 
when  a  small  fraction  of  another  gas  is  present.      A  possible  explana- 
tion on  the  basis  of  charge  exchange  and  the  metastable  states  of  helium 

[  1 8 1 

is  given  by  Sharpe  with  reference  to  experimental  findings 

NEGATIVE-ION  FORMATION 

An  important  gas  defect  is  negative-ion  formation  due  to  electron 
attachment  to  a  gas  molecule.  It  does  not  occur  in  pure  noble  gases 
nor  in  N?,  H  ,  CO,  CO?,  or  CH..  However,  it  is  likely  to  occur  in 
the  collisions  of  electrons  with  the  halogens,  0?y  or  water  vapor.  In 
the  development  of  Eq.  (Ill-  1 )  the  following  equation  was  derived  for 
the  case  where  a  fraction  h!  of  ionization  electrons  reaches  the  wire 
without  forming  negative  ions: 


i 

Vp 


(htmAm()) 


2    r 


V 

0     +  J_ 

2  mn 

mQ  0 


1.68 


(III- 2) 


where  V~  =  variance  of  number  of  ions  initially  liberated.      For  h1  -   1 
(no  negative-ion  formation),    this  equation  leads  to  the  previous 
approximation  that  relative  variance  is  about   l/mn  when  we  use  an 
experimental  result  that  V~       m_/3.      That  is,    the  relative  variance  is 
approximately 

\1 


J_ 
m 


0 


i 


/J 


(A  bibliographical  review  on  the  development  of  Eqs.  (Ill-  1 )  and  (III- 2 ) 
is  given  in  Reference  26). 

Thus,    for  accuracy  in  low-energy  work  it  is  recommended  that 
the  gas  be  purified  to  remove  negative-ion-forming  impurities.      How- 
ever,   since  the  factor  in  brackets  in  the  expression  for  relative 
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variance  is  multiplied  by  the  reciprocal  of  m        we  expect  that  in  our 

case  the  purification  requirements  may  not  be  so  stringent.      This  can 

r  ig] 
be  roughly  checked  by  an  approximate  expression  for  the  con- 

centration C  of  electronegative  gas  in  a  mixture  in  terms  of  the  gas 
pressure  P  (in  atmospheres),    and  Shubweg  6  (in  cm)  based  on  gas- 
discharge  theory: 

7  •     10 

C^-^r (m"3) 

The  Shubweg  6  is  the  product  of  the  mean  free  life  of  the  electron  and 
its  drift  velocity  under  the  applied  field.      Values  repre  seni  tative  of 
the  case  in  many  gases  were  assumed  for  the  electron  mobility  and 
probability  of  attachment  per  collision,    h,    for  the  gas. 

Suppose  we  wanted  a  gas  purity  that  would  make  h'  in  our 
expression  for  relative  variance  equal  0.99.      If  an  electron  started  a 
distance  r  from  the  central  wire  of  a  counter,    the  probability  of  its 

reaching  the  wire  without  attachment,    by  the  definition  of  the 

-r/6  -r/6  , 

Shubweg  6,    would  be  e  .      Therefore  we   set  e  =  0.99  or 

6     ~  100  r.      For  a  counter    in  which  3  cm  would  be  a  representative 

value  for  r  and  at  a  pressure  of  2.5  atoms  we  get,    from  Eq.    (Ill- 3), 

C   =  0.01%,    a  fairly  stringent  requirement.      If  we  try      h'  ~  0.9  and 

make  the  same  assumptions,    we  get  C     ~  0.1%,    or  99.9%  purity 

required.      Now,    going  back  to  the  approximate  expression  for  relative 

4 
variance  and  assuming  m~  ~   10    ,    we  can  solve  for  a  value  of  h*  that 

would  give  a  relative  variance  of,    say,    0.1%.      This  turns  out  to  be    a 

seemingly  very  easy  criterion,    hs  =  0.179.      However,    making  the 

same  assumptions  as  before,    we  find  C  is  about   1%.      But  in  the 

equation  for  C  the  assumptions  for  electron  mobility  and  h,    which 

exhibits  resonances  for  some  gases  dependtent  on  electron  energy, 

could  both  conceivably  be  off  by  almost  a  factor  of  ten,    and  considering 

also  other  approximations,    we  could  easily  be  back  to  our  99.99% 

purity  requirement.      The  parameters  in  the  gas-discharge    theory  are 

not  accurately  known  for  all  the  pure  gases  and  cannot  be  predicted 

for  mixtures,    especially  if  the  percentages  are  not  very  accurately 

known. 
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RADIOACTIVE  BACKGROUND 

The  problem  is  different  if  there  is  an  appreciable  tritium 
impurity  in  the  helium- 3  supply,      The  loss  of  energy  resolution  is  then 
due  to  the  pile-up  of  low- energy  background  counts  due  to  the  radio- 
active decay  of  the  tritium.      This  decay  has  an  end-point  beta  energy 
of  about   18.9  kev  L       J. 

To  estimate  the  required  purity  of  helium- 3  with  respect  to 
tritium  we  must  first  analyze  the  effect  of  pile-up  in  general 
Consider  the  piling  up  of  square  pulses  of  uniform  height  W  and  equal 
width    t  .      A  pulse  height  n  W  is  produced  if  n  pulses  occur  within  a 
time  t.      If  a  scaler  discriminatpr  level  were  set  at  n  Ws    an  erroneous 
count  would  then  be  registered.      We  can  get  approximate  values  for  a 
counter  by  taking  t  equal  to  the  equipment  resolving  time. 

If  the  occurrence  rate  of  events  in  the  counter  each  giving  rise 
to  an  energy  E  has  a  mean  rate  of  nn,    the  mean  number  of  events 
during  a  resolving  time  is  n^T  and  the  standard  deviation  is     ^nr«T  ^or 
n~T  »  1.      Thus,    the  energy  dissipated  per  resolving  time  is 
n^r  E±  n  n~  tE         Now  suppose  the  threshold  energy  of  the   system  is 
such  that  it  takes  an  energy  of  nrE  to  produce  a  spurious  count.      Then 
the  following  expression  gives  the  number  of  standard  deviations  by 
which  the  number  of  counts  during  a  resolution  time  must  exceed  the 
mean  number  in  order  for  a  spurious  count  to  be  produced: 

(nr  -   nQT)/  \l  nQT      -        [(n  -   nQ)  /  \/ n~  ]   n/"t~.  (III-4) 

From  this  we  can  estimate  the  probability  p     of  receiving  a  spurious 
count  during  the  resolving  time  for  a  given  n0  and  t.      The   spurious 
counting  rate  is  given  by  I     "■  p   /t.      If  we  estimate  r  as  0.5    jjisec  and 

s  & 

wish  to  keep  I     less  than  0.1  count  per  minute,    p     must  be  less  than 

-9  S 

I   t  =  5/6  •     10      .      This  probability  corresponds  to  a  deviation  above 

s 

the  mean  between  5  and  6   standard  deviations.      Let  us  require  that 
nn  be  less  than  n  by  6  standard  deviations.      The  average  energy  of  the 
tritium    beta  spectrum  is  about  5.7  kev.    therefore  let  E  equal  this 
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value.      The  range  of  beta  particles  of  this  energy  is   such  that  most 

of  the  beta  energy  will  be   spent  in  the  counter  gas,    as  shown  in 

3 
Chapter  VI.      Since  the  Q  of  the  He      (n,  p)  T  reaction  is  770  kev,    a 

value  of  nrE  less  than  this  will  not  be  detected  as  a  spurious  count. 

Of  course,    a  value  of  nrE  that  is  an  appreciable  fraction  of  770  kev 

destroys  the  accuracy  of  neutron-energy  determination  if  enough  of 

the  pile-ups  occur  in  coincidence  with  valid  neutron  counts.      However, 

the  number  of  accidental  coincidences  between  two  counting  rates  is 

proportional  to  the  product  of  the  rates     and  the  resolving  time;  in 

this  case  the  true  and  spurious  rates  are  expected  to  be  small  enough 

that  the  problem  of  modification  of  true  pulse-height  counts  will  not 

be  as  serious  as  the  number  of  spurious  counts  recorded  compared  to 


the  number  of  true  counts.      If  we  set  iitE 
v8 


10v 


Using  Eq.    (III-4),    we  obtain 

3 


108  -  n 


0 


n 


0 


ITT 


(0,5 


770  kev,    n  is  approximately 

io-V/z. 


,8 


Solving,    we  have  n~  E   1.86  •     10     disintegrations  per   second. 

If  this  is  the  allowable  mean  disintegration  rate,    we  can  solve 
for  the  allowable  number  of  tritium  atoms  in  the  counter  from  the 
disintegration  equation: 


dN/dt  =     -    (0.693/T,)  N  -     -  nQ 


where 


N 
T 


This  gives  N  =   1.03 


1/2 
10 


=    number  of  tritium  atoms 

••  tritium  half  life  =   12.Z6  years 
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tritium  atoms.     Or  if  each  molecule  is 
1,6 


diatomic.,    the  allowable  number  is   5.65  -    10        molecules.      The 

sensitive  volume  of  the  central  counter  is  V     "  888  cc,    and  at  P  =  Z.5 

C  22 

atmospheres,    the  total  number  of-molecules  is     D  P  V     =  5.66  ■     10 

o  C  19 

molecules,    where  D  -  Loschmidt's  number  at   15     C,    or  2.55  '    10 

This  gives  a  tritium  composition  of  about  0.0001%.      This  looks  like  a 
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very  serious  purification  problem  if  there  is  significant  tritium  contam- 
ination in  the  helium-3  supply. 

SATURATION  EFFECTS 

It  is  apparent  that  for  high-energy  particles  and  high  enough 
multiplication,    saturation  effects  near  the  wire  should  eventually 
destroy  the  proportionality  of  the  multiplication  process.      Experiment- 
ally,   the  ratio  of  pulse   sizes  has  been  found  accurate  to   1.5%  for  M 

[28l 
less  than  a  critical  value  of  M    (E)  at  which  saturation  effects  begin 

For  the  counter  used  the  experimenters  found 

E  •    M    (E)  =  constant  =   108  ev.  (III-5) 

Using  this  for  our  counter  application,    where  the  maximum  (E     +  Q)  is 
1.77  Mev,    and  allowing  a  safety  margin  of  a  factor  of  two,    we  have 

M     (E  =   1.77  Mev)     *    28 
c 

If  we  also  check  expected  charge-to-length  ratio  in  avalanche  compared 

to  charge-to-length  ratio  on  the  wire,    a  ratio  of  4%  corresponding  to 

Eq      (III- 5),    it  appears  that  adoption  of  the  equation  for  this  counter 

should  ensure  no  significant  loss  of  proportionality  due  to  saturation. 

High  values  of  multiplication  are  also  not  desired,    owing  to  increased 

voltage- stability  requirements,    because  the  curve  of  multiplication 

versus  counter  voltage  in  most  cases   starts  out  rather  flat  but  soon 

rises  into  a  steep  ascent. 

POSITION  OF  IONIZATION 

Because  of  the  long  travel  time  of  positive  ions  to  the  cathode, 
amplifier  time  constants  must  be  used  that   ignore  part  of  their  con- 
tribution to  the  pulse.      This  makes  the   size  of  the  pulse  dependent  on 
the  position  of  ionization.      A  general  rule  is  to  use  M    >  10  to  avoid 
this  effect  when  counter  radius  is   100  times  wire  radius.      Let  us 
investigate  the  effect  for  our  counter,    for  y  =  b/a  =   1000,    and  for  no 
multiplication.      On  the  basis  of  the  expressions  for  fraction  of  total 
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induced  voltage  due  to  electrons  and  positive  ions,    Staub  shows  that 
the  fraction  of  the  volume  of  the  chamber  in  which  the  ionization  can 
originate  so  that  more  than  half  of  the  pulse  is  caused  by  the  fast 
electron  collection  is  b/(b+a),    since  the  contribution  to  the  pulse  is 


evenly  divided  when  r  =     M  (a  •    b)    is  the  position  of  original  ioniza- 
tion *■         .      Or  this  is  y/^y-fl),    which  equals   1000/1001,    which  seems  to 
indicate  that  multiplication  might  not  be  necessary  for  large  y  values. 
This  assumed  position  corresponds  to  r  s  if. 6a  for  our  counter.      As 
a  further  check  we  assumed  the  ionization  took  place  at  r  =  316  a, 
and  calculated  the  fraction  of  the  pulse  due  to  electron  motion 
(ionization  in  about  90%  of  the  volume  would  give  larger  pulses), 
using  the  assumptions  about  ionic  motion  used  in  Chapter  II  and 
some    other  equations  from  Reference   18.      This  fraction  was  0.835. 
Amplifier  time  constants  in  reality  would  be  long  enough  to  register 
some  of  the  positive-ion  motion,    which  would  improve  this  fraction. 
However;    starting  from  r  =  316a  under  the  same  conditions  and 
assumptions,    we  calculate  the  ionic  motion  to  contribute  only  about 
0.00Z  to  the  total  pulse  in  a  collection  time  of  5  fj.sec.     Thus  we  con- 
clude  some  multiplication  is  necessary  to  remove  effects  due  to 
position  of  ionization. 

END  EFFECTS 

With  no  preventive  measures  there  would  be  an  abnormally 
high  electric  field  at  the  ends  of  a  counter  wire.      The  loss  in 
resolution  thereby  can  be  partially  erased  by  having  the  wire 
shielded  or  thickened  near  the  ends  to  give  a  low  field.      The  variation 
of  M  along  a  counter  wire  has  been  measured  and  found  to  be  import- 
ant at  a  distance  on  the  order  of  a  counter  radius  from  the  thickened 

[291 
ends  .      In  a  high-pressure  proportional  counter  designed  to  be 

used  for  neutron  spectroscopy  it  was  found  that  the  pulse-height 
distribution  varied  radically  over  the  extent  of  the  counter  in  a  pre- 
liminary design  without  field  tubes'-         .      The  loss  of  resolution  com- 
pared with  the  resolution  at  the  center,    was  significant  when  the 
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investigators  moved  a  narrow  pencil  of  thermal  neutrons  with  which 
they  were  probing  the  counter  0.5  inch  away  from  the  center.      The 
wire  was  about  8  inches  long.      The  simplest  way  of  overcoming  end 
effects  is  to  use  a  long  counter- -that  is,    long  compared  with  the 
radius.      Of  course,    this  is  not  always  practicable.      Our  counter, 
with  wire  length  21.6  cm  and  counter  diameter   7.62  cm,    is  not  a 
very  long  counter,    and  we  might  expect  considerable  loss  of  resolu- 
tion due  to  end  effects.      The  end  effects  can  be  eliminated  by  mechan- 
ical methods,    which  add  complexity  to  the  counter.      A  divided  counter 
was  constructed  in  which  a  bead  divided  the  wire  into  two  unequal 
lengths1         .      Subtracting  one  spectrum  from  the  other  eliminates 
end  effects.      Field  tubes  extending  over  a  guard  tube  eliminate     end 
effects  when  field-tube  potential  is  adjusted  to  a  proper  value  and 

field  and  guard  tubes  are  correctly  located.      Field  tubes  were  proposed 

f  3  1  "I  f32l 

in   1951  and  used  the   same  year1  to  investigate  the  low-energy 

spectrum  of  tritium  down  to  approximately  200  ev. 

WALL  EFFECTS 

Wall  effects  as  outlined  for  a  gas  recoil  chamber  by  Rossi  and 

[21] 
Staub  are  pertinent  for  charged-particle  reactions,    since  in  both 

cases  the  loss  in  resolution  results  when  the  particle  that  causes  the 

ionization  in  the  gas  does  not  expend  its  full  range  in  the  sensitive 

volume.      Three  cases  are:     (a)    The  reaction  takes  place  in  the 

sensitive  volume  but  one  of  the  charged  particles  produced  crosses 

the  end  boundary;     (b)     The  reaction  takes  place  in  the  sensitive 

volume  but  the  charged  particle  hits  the  lateral  wall;     (c)    The 

reaction  takes  place  outside  the   sensitive  volume,    but  one  of  the 

charged  particles  enters  the   sensitive  volume.      The  anticoincidence 

ring  does  not  help  in  those  instances  in  which  the  charged  particle 

crosses  the  ends.      Thus,    it  is  desirable  to  have  a  small  ratio  of  area 

of  ends  to  lateral  area,    or-- what  is  the   same  thing- -a.long  counter. 

as  in  reducing  end  effect.      Pulses  reduced  in  size  by  the  wall  effect 

contribute  a  continuum  of  energies  from  zero  to  (E     +  Q)  to  the 
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pulse-height  distribution.      Without  the  use  of  an  anticoincidence  ring  or 
a  magnetic  field  to  curve  the  paths  of  the  charged  particles,    the  wall 
effect  can  be  reduced  only  by  increased  size  of  the  counter  and  in- 
creased gas  pressure. 

OTHER  EFFECTS 

Inaccuracies  in  energy  resolution  can  also  come  from  variation 
in  applied  voltage,    variation  in  position  and  diameter  of  the  wire,    and 
dust  particles  on  the  wire.      These  are  discussed  in  the  general 
references.      Any  alone  could  be  very  detrimental,    but  it  appears 
that  all  can  be  avoided  with  careful  techniques. 
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CHAPTER  IV 
GAS  FILLINGS 

The  counter  was  first  tested  by  filling  from  a  bottle  of   argon 
plus  4%  carbon  dioxide.      The  filling  and  evacuation  system  used  had 
an  oil  diffusion  pump  and  liquid  nitrogen  trap,    but  no  special  purifica- 
tion procedure  or  trap  was  used  in  the  gas-filling  line.      Satisfactory 
operation  was  obtained  in  counting  gamma  rays  at  pressures  of  one, 
three,    and  five  atmospheres.      Since  the  secondary  electrons  have  a 
continuous  distribution  of  pulse  heights,    we  would  not  expect  a  flat 
plateau  in  the  curve  of  counting  rate  versus  voltage  in  the  proportional 
region,    but  we  do  desire  a  slope  that  does  not  place  too  exacting 
requirements  on  voltage  regulation.      At  5  atmos  we  found  about    -|% 
increase  in  counting  rate  per  volt  in  the  2600-to-  2700-  volt  range  on 
the  main  counter.      This  was  at  a  high  value  of  M;  the  range  of  pulse 
heights  available  was  several  thousand  times  the  minimum  detectable 
pulse  height.      At  3  atmos  in  the   same  voltage  range,    the  change  in 
counting  rate  per  volt  was  about  3%.      Of  course,    the  M  was  greater 
owing  to  the  lower  pressure.      The  ring  counter  at  5  atmos  had  about  1% 
increase  per  volt  in  the  2  100-to- 2200-volt  range.      Its  counting  rate 
was  higher  owing  to  its  proximity  to  the  wall. 

The  helium  filling  was  from  a  bottle  of  99-99%  pure  helium;  the 
desired  percentage  of  carbon  dioxide  was  previously  introduced  through 
a  separate  filling  line.      Both  filling  lines  were  run  through  dry  ice- 
acetone   slush  traps  when  CO,  was  used  or  a  liquid  nitrogen  trap  if 
helium  alone  was  used.      (CO,  would  be  frozen  out  by  a  liquid  nitrogen 
trap.  )     The  use  of  the  traps  did  not  seem  to  have  a  large  effect.      The 
CO,  was  found  very  necessary  to  avoid  spurious  discharges  in  the 
counter . 

The  stabilizing  effect  of  polyatomic  gases  mixed  with  noble  gases 
seems  to  depend  on  their  ability  to  absorb  photons,    in  contrast  to  the 
small-photon  capture  cross   section  of  the  noble  gases   (many  of  these 
photons  are  produced  by  metastable  atoms),    and  also  possibly  on  their 


25 


ability  to  receive  energy  from  metastable   states  that  have  high 

excitation  energies       Colli  has  investigated  ultraviolet  photons  in  the 

decay  of  metastable  argon  atoms  and  found  them  to  be  responsible, 

through  a  photocathodic  process,    for  the   starting  and  maintenance  of 

["33] 
corona  current  in  cylindrical  argon  counters  .      The  photons  were 

produced  in  the  Townsend  avalanche  on  the  wire.      Another  effect  of 

metastable  atoms  is  the  release  of  secondary  electrons  from  metal 

[34] 
surfaces  .      The  positive  ions  can  also  eject  secondary  electrons 

from  metals,    the  requirement  being  that  the  energy  of  ionization  be 

greater  than  twice  the  work  function  of  the  metal,    which  is  about  5  ev 

for  common  cathode  walls.      Helium  is  expected  to  be  more  troublesome 

than  other  gases  because  of  its  high  ionization  potential,    about  24.5  ev, 

and  the  high  energies  of  its  metastable  atoms,    about  20  ev.      Evidence 

for  the  difficulty  of  interpreting  experiments  with  helium  has  been 

[35] 

presented  on  the  basis  of  the  existence  of  the  helium  molecule-ion 

and  metastable  molecule       It  had  been  shown  the  He       ion  is  present 
at  very  low  pressures,    and  the  concentration  increases  with  pressure 
relative  to  the  concentration  of  He    . 

The   stabilizing  effect  of  CO?  was  demonstrated  in  the  counter  by 
determining  the  approximate  maximum  voltage  on  the  main  counter 
that  could  be  applied  without  causing  a  rapid  breakdown  when  a  radio- 
active source  was  placed  a  given  distance  from  the  counter.      Total 
pressure  was   5  atmos.      The  results  are  listed  in  Table  II.     The  break- 
down   caused  by  background  radiation  alone  occurs  with  pure  helium 
at  800  volts. 

However,    adding  more  CO-,  was  not  a  solution  because  the 
pulse  heights  were  thereby  drastically  reduced.      It  had  been  found 
previously  that  in  mixture  of  argon  plus    10%  CO?  there  was  a  loss  of 
pulse  height  above  pressures  of  3  atmos,    probably  due  to  columnar 
recombination,    which  is  high  in  CO?.      This  refers  to  recombination 
along  the  track  of  a  densely  ionizing  particle.      It  was  not  possible  to 
attain  a  large  enough  range  of  pulse  heights  above  a  minimum  dis- 
criminator  setting  for  counting  recoils  from  a   Pu-Be   (a,  n)  neutron 
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Table  II.      Percentage  of  CO,  in  helium  versus  voltage 
for  breakdown  on  main  counter. 

Percentage  of  CO  Critical  Voltage 

4  1500 

6  2000 

12  2500 

20  2800 

25  3000 

33-1/3  

source  without  going  to  a  counter  voltage  at  which  spurious  discharges 
would  occur.      The   situation  improved  with  decreasing      pressure. 

A  final  filling  at  2.5  atmos  with  2%  CO     was  chosen  for  further 
experimentation       The  curve  of  counting  rate  versus  voltage  for  count- 
ing neutron  recoils  from  the  Pu-Be  source   showed  a  slope  of  0.84% 
change  in  counting  rate  per  volt  in  the  range  of  1100  to   1200  volts  on  the 
main  counter  wire.      For  the  outer  wire  there  was  a  4.3%  change  in 
counting  rate  per  volt  in  the  range  from  800  to   1000    volts.      The 
voltage  -  stability  requirements  at  5  atmos  at  a  voltage  at  which  the 
spurious  discharges  did  not  occur  were  not  significantly  different. 

One  advantage  of  helium  as  the  counting  gas  is  readily  apparent- - 
its  reduced  sensitivity  to  background  gamma  radiation.      The  range  of 
the   secondary  electrons  in  g/cm     is  about  the   same,    but  helium  has 
about   l/lO  the  density  of  argon.      The  low  atomic  number  of  helium  is 
also  an  advantage  in  that  the  interaction  of  gamma  rays   directly  with 
the  gas  is  small. 
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CHAPTER    V 
COUNTING  RATE  AND  EFFICIENCY 


ing  rate  may  be  expressed         as 


For  an  isotropic  neutron  flux  in  a  proportional  counter  the  count- 
expressed         as 

I-NVc    [   a(En)    c(>E     F(En)     f(En)dEn, 

where 

N  =  density  of  reacting  nuclei  (nuclei /cm    ), 

3 
V     '  volume  of  absorbing  material  (cm    ), 

<j  (E    )  -  neutron  cross  section  of  absorbing  nuclei 
n  & 

(cm    ), 

4> -rp      =  neutron  flux  per  energy  interval  (neutrons/cm    - 
sec-erg), 

F(E    )  =  detection  efficiency  defined  as  ratio  of  count- 
ing rate  to  rate  of  reacting  events. 
Ideally,    for  an  exo$^rgic  reaction  and  integral  detection,    F(E    )  is  a 

step  function:  zero  for  E      <  B  and  unity  for  E       >    B,    where  B  is 
r  n  7  n 

threshold  bias  energy  for  equipment.      This   step  function  is  modified 
by  a  geometric  factor,    G(E    ),    for  wall  effect  as  shown  below.      For 
differential  detection.,    F(E    )  is  ideally  a  difference  between  two  step 
functions  with  thresholds  B,   and  B?.      For  a  recoil-detection  process 
in  which,    with  no  collimation,    a  neutron  of  any  energy  may  give 
pulses  from  zero  to  a  maximum,    we  must  evaluate  F(B/E    )s    which 
gives  the  fractional  number  of  recoils  giving  pulses  larger  than  B. 
Thus,    for  a  recoil  detector,    one  has 

F(E    )  fc  G(E    )  F(B/E    ). 

n'  v    n'       v    '      n! 

The  term  f(E    )  above  is  a  correction  factor  for  perturbation  of 
n'  r 

flux  caused  by  insertion  of  the  absorber.      This  factor  is  neglected  in 
the  rest  of  the  discussion. 

A  solution  of  the  equation  for  counting  rate  would  in  most  cases 
be  extremely  complicated.      Some  insight  can  be  gained  by  considering 

the  simpler  case  of  a  monoergic  neutron  flux.      Then  we  can  eliminate 

2 
the  differential  and  replace  <f>F     with  <j>  (E    )  in  neutrons/cm    -sec. 

n 
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For  neutrons  of  a  given  energy    this  leads  to 

I  =  N  V     ct(E    )    <b(E    )     F  (E    )  . 
c  n      T      n  nf 

The  efficiency  of  a  counter  is  defined  as  i/q,    where  q   is  the 

number  of  neutrons  per  second  traversing  the  counter  volume; 

q  =  cj)(E    )A,    where  A  is  the  effective  exposed  area  of  the  counter.      Then 

counter  efficiency  is 

NV     a(E    )    <b(E    )  F(E    )  NV     a  (E    )  <b(E    ) 

c  n      Y  v    n'       x    n  c  n    Y      n' 

e. 


C  (b(E    )  A  A 

If  we  call  V  /A  =  d,    the  average  distance  traveled  by  a  neutron  in 

•        .v,C  «.      [24]  ,  „.    • 

crossing  the  counter  ,    counter  efficiency  is 

er   =  N  a    (E    )  d  F(E    ]  =  e0   F(E    )  , 
C  n'  v    n'  R      v    n'  ' 

where      e_   -  N  a    (E    )  d  is  the  efficiency  of  radiation  defined  as  the 
R  n'  y 

ratio  of  the  rate  of  reacting  events  to  rate  of  neutrons  traversing  the 
counter  volume. 

The  quantity  N  is  readily  calculated  by  expressing  it  as   D  P, 
where  D  is   Loschmidt's  Number  and  P  is  the  pressure  in  atmos. 

For  the  counter,    we  have  d  =   1,59  cm  if  we  take  the  total  sur- 
face area  of  the  main  counter  as  the  effective  exposed  area  for  an 
isotropic  flux.      If  we  assume   Z.5  atmos  of  helium-3  and  a  neutron 
energy  for  which  the  (n,  p)  cross  section  is  one  barn,    we  have 

e_   ~    0.01%.      For  the  counter  filled  with  ordinary  helium  at  the  same 
is. 

pressure  and  for  a  neutron  energy  for  which  the  scattering  cross 
section  is  2  barns,    we  have  eR  ~  0.02%  for  recoils. 

The  opposite  extreme  to  an  isotropic  flux  is  a  collimated  neutron 
beam.      Here  the  general  equation  for  counting  rate  is 


I  =  A  (  <>E       (1   -   e"N    a(En}  T)  F(En)  dE{ 
n 


where 


A  =  area  perpendicular  to  beam, 

T  =  thickness  in  cm  of  material  parallel  to  beam. 
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Assuming  a  monoargic  neutron  beam  as  before,    we  can  write 
I  =  A^En)[,.e-N-(E„»T]     F(En), 

ec,[l-e-N-<En»T]F«En,  =  eRr(En)> 

and  for  N  a   (E    )  T    «     1,    e„  «    N    c   (E    )  T. 
n'  R  v     n' 

In  this  case  we  see  that  the  radiation  efficiency  of  a  cylindrical  counter 

in  a  neutron  beam  is  dependent  on  its  orientation  with  respect  to  the 

beam.      If  the  beam  direction  is  perpendicular  to  the  long  axis  of  a 

counter  of  radius  b,    solving  for  the  average  thickness  of  path  gives 

(tr/Zlb  or  5.98  cm  for  the  main  counter.      Thus,    the  counter  has  a 

higher  radiation  efficiency  under  these  conditions  by  the  factor 

T/d  =  3.76.      Of  course,    e_   would  be  even  higher  for  a  beam  down 

the  long  axis  of  a  cylindrical  counter. 

EXPERIMENTAL  SETUP 

However,    in  experimental  setups  the  counter  usually  is  neither 
in  an  isotropic  flux  nor  exactly  in  a  beam.      Consider  placing  a  neutron 
source  at  point  P,    a  distance  c  from  the  counter  axis  and  opposite  the 
center  of  its  length,    as   shown  in  Fig.    4.      If  c  is  not  too  small  and  the 
counter  dimensions  not  too  large,    it  seems  we  could  approximate  the 
expected  counting  rate  by  assigning  an  average  neutron  flux  to  the 
counting  volume: 

4>        -  Qe/4irr2    ~  , 

av  s'  elf 


where  Q     =  rate  of  emission  of  neutrons  by  the   source  and 

J f~2 2 2 

r    ,,  =   si  (c  d)  ,    where  d  =  sib     +  f     +  c      .      Then, 
eff 

I  =  N  V     a  Q/  4tic  d 
c  ' 

(assuming  F(E    )  =   1  for  the  time  being).      The  rate  for  neutrons 
traversing  the  counter  in  neutrons  per  secrond  can  then  be  expressed 
as  this  average  flux  times  the  area  of  the  central  rectangle,    or 
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Fig.    4.      Geometrical  arrangement  of  neutron  source  and 
counter. 
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.zb.z£=    Q*         « 


4-rrcd  it         •      cd 

This  expression  was  set  down  to  show  that  this  method  is  equivalent 

to  a  common  method  of  computing  the   solid  angle   subtended  by  a  finite 

r  i  o  i 
rectangular  counter  .      For  one  of  the  small  rectangles  of  area  b  f 

in  the  geometry  shown  the  equation  for  solid  angle  given  is 

tan  ft  =  b  f/c  d 
For  a  small  angle  we  have  tan  ft~ft  and  for  the  entire  rectangle  the   solid 
angle  would  be 

4  b  f/c  d. 
Thus,    we  have 

q  Tc7    /    4T    Qs"(Qs/tt)'    <bf/cd)  , 
as  above.      Also  one   should  note  that  we  have  assigned  the   same 
radiation  efficiency  in  this  geometry  as  for  a  beam,    as  we  can  show 

e       =  i/q  =  N  a   Vc/4  bf  (still  neglecting  F(E    ). 
But  we  have  V  /4  b  f  =  ub^  •    2  f/4  b  f  =  (TT/2)b,    the  average  thickness 
found  for  a  beam  through  the  counter. 

An  arrangement  like  this  was  used  with  a  Pu-Be  neutron  source 
of  Q     =   1.54  •     10     neutrons/sec  and  d  =  21.2  cm.      Then  we  compute 
l/F(E    )  =   1662  cpm  for  a  pressure  of  2.5  atmos  of  helium,    and  using 
an  average   scattering  cross  section  of  2  barns  corresponding  to  an 
average  neutron  energy  of  about  4  Mev.      (Actually  the  Pu-Be  neutron 
spectrum  is  complex  from  zero  to  about   11  Mev,    and  no  attempt  has 
been  made  to  compute  an  accurate  average  value.  ) 

If  we  obtain  curves  at  different  values  of  high  voltage,    as   shown 
in  Fig.    5,    we  can  estimate  the  point  where  they  tend  to  zero  pulse 
height.      At  this  point  F(E    )  is  unity  and  the  experimental  value  of 
counting  rate   should  correspond  to    1662  cpm.      It  appears  to  be  only 
about  400  cpm,    not  a  close  correspondence,    but  there  were  so  many 
inexactly  known  quantities  and  dubious  averaging  processes  that  it  is 
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Fig.    5.      Integral  bias  curves  for  main  counter  at  various 
counter  voltages, 

A.  1200  volts, 

B.  1165  volts, 

C.  1135  volts, 

and  using  neutron  source  in  arrangement  of  Fig.    4. 
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hard  to  draw  definite  conclusions.      The   scattering  cross  section  is 
rapidly  varying  in  this  energy  range  and  the  manner  in  which  the  inte- 
gral bias  curve  approaches  zero  is  not  known.      Also  consider  the  count- 
ing rate  of  175  cpm  for  the  highest  counter  voltage  used  at  fc.  pulse 
height  of  10  volts,    the  lowest  discriminator   setting  for  which  counts 
were  recorded.      If  we  let  G(E    )  =  0.83--a  value  computed  for  an 
average  alpha  recoil  in  Chapter  VI  then,    for  a   10-volt  biass    we  have 

F(B/E    )  =  0.127,    which  would  indicate  a  detection  of  a  small  fraction 
'      n 

of  the  recoil  alpha  particles.      The  curves  for  the  ring  counter,    Fig.    6, 
seem  to  tend  to  a  much  higher  zero-point  counting  rate,    but  it  is 
believed  that  other  events  than  the  recoils  are  being  counted,    since  the 
volume  of  the  ring  counter  is  less  than  that  of  the  central  counter. 

RECOIL-SPECTRUM  CONSIDERATIONS 

It  should  be  much  easier  to  draw  some  definite  conclusions  about 
the  suitability  of  the  tube's  multiplication  factor  for  obtaining  usable 
pulse  heights  if  a  monoergic  neutron  source  were  used,    preferably  one 
giving  a  recoil  alpha  in  the  energy  range  of  interest,    770  kev  to   1.80 
Mev.      The  counting  results  from  a  monoergic   source  would  be  much 
easier  to  analyze. 

This  is  illustrated  as  follows.      For  a  monoergic  source  the 

neutron  spectrum,    N(E    )  versus  E    ,    is  ideally  a  spike  and  N(E    )  is 

the  total  number  of  neutrons;  the  differential  pulse-height  spectrum 

due  to  the  recoils,    n(W)  versus  W,    where  W  is  pulse  height,    is  a 

rectangle;  and  the  integral  bias  recoil  spectrum,    N(W)  versus  W,    is 

ideally  a  constant- slope  line.      These  are   sketched  in  Fig.    7.      These 

shapes  of  the  differential  and  integral  bias  curves  are  justified  as 

follows.      The  pulse  height  depends  on  the  energy  of  the  recoil  nucleus; 

the  counter  and  amplifier  convert  that  energy  into  pulse  height.      The 

recoil  nucleus,    as   shown  in  Chapter  VIII,    can  have  a  maximum  energy 

of  K     E    ,    where  K     -  4  A/(l   +  A)     ,  if  we  usel  mass  numbers  for  the 
an  a  ' 

masses.      Let  W  be  the  energy  of  recoil  nucleus   (converted  into  volts 
of  pulse  height  by  the  counter  system).      Then  W  can  vary  from  0  to 
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Fig.    6.     Integral  bias  curves  for  ring  counter  at  various 
counter  voltages, 

A.  1025  volts, 

B.  975  volts, 

C.  925  volts, 

and  using  neutron  source  in  arrangement  of  Fig.    4. 
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Fig.    7.      Ideal  curves  for  monoenergetic  neutron  spectrum. 

A.  Neutron  spectrum. 

B.  Differential  bias  recoil  spectrum. 

C.  Integral  bias  recoil  spectrum. 
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K     E     and  the  number  of  pulse  heights  from  recoils  from  the  monoergic 

source  are   spread  out  over  this  range  of  W  values.      The  question  is 

whether  or  not  the  distribution  is  uniform  over  the  range,    giving  a 

rectangle.      At  the  values-  of  E     of  interest  here,    and  for  light  nuclei, 

it  is  accurate  to  assume  that  neutron  scattering  is  mainly  S-wave  and 

spherically  symmetric  in  the  center- of-mass   system  of  coordinates. 

We  wish  then  to  express  the  recoil  energy  in  terms  of  E     and  the 

angle  of  scattering  in  center- of-mass  coordinates   (both  energies  in 

laboratory- system  coordinates).      Using  conservation  of  energy  and 

r  o/  "I 
an  equation  developed^  for  relating  the  velocities  of  the  neutrons 

in  lab  coordinates  before  and  after   scattering  with    9,    the  angle  of 

scattering  in  cm.    coordinates,    we  can  show  W  =  \.YL   E      (1-cos  Q): 

°  2     a    n 

then  dW  =   4  K     E     sin  0  dO.      But  in  spherical  coordinates,    we  have 

2      a      n  r 

dB  =  dft/2-rr   sin    0  ,    where     dQ,  ~  element  of  solid  angle,      therefore 

dW  =  K     En  (dn/4'ir).      We  can  say      Prob(W)  dW  =  (a  (0)/ct    )  dfi  , 

where 

g(6)  =  differential  scattering  cross  section  (cm.  ), 

<7  -  integral  scattering  cross  section, 

Prob(  W)  dW  =  probability  for  recoil  nucleus  to  acquire  an  energy 

between  W  and  W  +  dW  , 

(cr(6)/(J    )dQ,  ~  probability  that  neutron  will  be  scattered  through 

angle     6     into  d£2  . 

It  follows:      Prob   (W)  -  a  (Q)/ a      -      4tt/K     E     .  For  an  isotropic 

v     '  v    "      s  'an  ^ 

scattering  distribution  (a  (  6)/ o      "•     1/4tt)  and  monoergic  neutrons, 
Prob(W)  is  thus  a  constant,     l/K     E    ,    and  n(W)  is  a  rectangle  of 
constant  value  proportional  to  N(E    ),      The  integral  bias  distribution 
is  an  integral  of  the  rectangle,    integrating  from  maximum  W  to  zero, 
giving  the  triangular  shape. 

For  a  non-monoergic  neutron  spectrum,    N(E    )--instead  of  being 
the  total  number  of  neutrons- -represents  the  number  of  neutrons  per 
energy  interval,    and  we  must  integrate  over  the  spectrum.      We  can 
say 
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n(W)  :-  Kj  [   Prob     W    N  <En\    (Js   (En)  dEn    -  K{  (  l/K^       N-E^  ^En)dEn 

where  K,   is  an  arbitrary  constant  and  a        E    1  is  the  elastic  scattering 
1  s        n  e 

cross  section  as  a  function  of  neutron  energv       Thus     for  a  complex 

neutron  spectrum  the  shapes  of  the  differential  and  integral  bias 

curves  are  very  difficult  to  predict.      For  performing      the  integration 

in  analytical  form  an  expression  for  NsE    )  and  a    (E    Imust  be  usfd 
7  r  n  s       n 

Generally,    for  most  sources  N 'E    }  goes  from  zero  at  E  0  to  a 

n"  °  n 

maximum  and  back  to  zero  at  E     -  E    (max).      The  intervening 

n  nx  6 

structure  may  be  quite  complicated.      Let  us  assume  a        E    )  varies  as 

j-  .  g  ]  s        n 

1/W  E  .      This  is  approximately  true1  for   scattering  on  protons 


for  E     between  0.2  and  5  Mev      an  approximate  relation  being   (.a    )     ~ 
0 —  s  p 

4.5/mE     for  E     in  Mev  and    a      in  barns.      Et  is  less  true  for  helium 
'  n  n  s 

which  has  a  peak  of  about  b.b  barns  at    1.15  Mev  and  then  tails  off 

Then  if.    for  illustrative  purposes,    we  assume 

N(E    )  -  E    3/Z     sin  E    , 
n"  n  n 

where  E     goes  from  zero  to  E      (max)  ~    tt     ;.n  arbitrary  units     we  have 
n  °  n 

a   spectrum  shape  that  can  be  readily  integrated: 

E      =TT 

n 

n(W)  =  K,         f         l/(K     E    )     E    3' 2  sm  E     E    "i/2dE 
v  1       J  'ann  nn  n 

E    =W/K 


K 


i 


K 
a 


(1  +  cos   (W/K  )), 


where   W  goes  from  0  to  K     tt.      To  save   writing   W/K     as  the  variable 

el  cL 

on  the  next  integration,    let  us  change   W  from  energy  to  arbitrary 
voltage  units  where   W  goes  from  0  to  tt,      Since  we  have  divided  the 
W  scale  by    K    .    we  must  multiply  the  density  d  s;r  bution  by  K    .      The 

3.  cL 

differential  distribution     becomes 

n(W)  =  K.      I  +  cos  W 
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The  integral  bias  distribution  is  obtained  by  integrating  this  from 
the  right:  w 

N(W)  =  -  K  [    (1  +  cos   W)  dW  =  Kx   (tt  -    W  -   sin  W). 

TT 

The  value  of  this  distribution  at  W  -  0  is  the  total  number  of  scattering 
events,    ttK,.      This   should  equal  the  number  of  neutrons  times  the 
probability  of  scattering  or 


TT  TT 

K.      (    N(E    \  <j      (E    )  dE     =  K.        f    E      sin  E     dE     =  ttK.  . 
llnsvn'n  ljn  nn  1 

0  0 

The  three  distributions  are  sketched  in  Fig.    8.      The  significant 
trend  is  the  way  the  integral  bias  distribution  rises  rapidly  near  zero 
pulse  height,    even  with  this  assumed  neutron  spectrum  weighted 
heavily  toward  the  high-energy  end  of  the  spectrum.      With  a  complex 

neutron-energy  spectrum--as  obtained,    for  example,    from  (a,  n) 

[371 

sources  --the  shape  of  the  integral  bias  curve  and  the  manner  in 

which  it  tends  to  the  zero  pulse-height  counting  rate  becomes  very 
difficult  to  predict. 
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Fig.    8.      Assumed  neutron  spectrum  and  recoil  spectrum 
for  assumed  conditions. 

A.  Neutron  spectrum. 

B.  Differential  bias  recoil  spectrum. 

C.  Integral  bias  recoil  spectrum. 
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CHAPTER    VI 

PARTICLE  RANGE  AND  WALL  EFFECT 

Knowing  the  range  of  an  alpha  particle    or  proton  in  air,    we  can 
compute  the  range  in  helium  if  we  know  the  relative  atomic  stopping 
power  of  helium  compared  with  air.      The  range-energy  curves  to  be 

used  for  alpha  particles  and  protons  in  air  at   15     C  and  760  mm 

[38l 
pressure  are  those  of  Bethe,  which,    on  the  basis  of  experiment, 

corrected  for  low  energies  the  original  curves  by  Livingston  and 

[39] 

Bethe  ,      The  original  work  was  a  theoretical  treatment  based  on  the 

Born  approximation  but  also  modified  by  experimental  data  then 
available. 

The  basic  equation  giving  energy  loss  per  cm  of  path  defines 
stopping  power: 

-  dE/dR  =     (4it  e4z2/mv2)  NB, 
where 

E  =  energy  of  incident  particle, 

R  =  path  length, 

e    =  electron  charge, 

z    =  atomic  number  of  incident  particle, 

m  =  electron  mass, 

v    =  velocity  of  incident  particle, 

N  =  number  of  atoms  per  cc  of  stopping  material, 

B  =  stopping  number  =  Z  In  (2m  v   /I      )a    where 

Z  =  atomic  number  of  stopping  material, 

I        =  K  Z  =  average  excitation  potential  of   atom  of  stopping 

material, 
K  =  Bloch  constant  (which  is  not  actually  a  constant  but    varies 

somewhat  for  different  atoms  and  is  usually  determined 
experimentally  because  it  is  difficult  to  calculate 
theoretically). 
The  relative  stopping  power,    compared  with  air,    is  defined: 
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s  =  B/BQ  -  Z/Z0[(ln  2mvZ  -  In  I      }/<ln  2mvZ  -  In  IQ)]    , 

where  the  zero  subscripts  refer  to  air.      Here   s  is  seen  to  be  a  function 

of  incident-particle  velocity     and  for  different  atoms  increases  more 

slowly  than  Z  alone.      We  find  it  accurate  enough  for  the  purpose  at 

hand  to  use  this  equation,    although  many  modifications  to  the  theory 

can  be  made  in  computing  B     as  done  by  Hirschfelder  and  Magee^ 

for  several  different  substances  for  protons. 

We  take  experimentally  determined  values  of  air  from 

Reference  39:  I«  -  80.5  ev  and  Z„  =  7.22.      For  helium,    I        =  44  ev 
U  0  av 

The  calculation  shows,    for   1-Mev  protons,    sTT     -    0,328.      To  convert 

r  He 

the  range  in  cm  in  air  at  NTP  (15     C,    760  mm)  to  that  in  helium  under 
the  same  conditions  we  note  that  stopping  power  is  proportional  to 
N  and  B,    therefore  the  range  in  cm  must  be  inversely  proportional  to 
these  quantities.      The  range  in  helium  is 

R  =  RQ  •    NQ/N  ■    s0/sHe  , 

where  the  zero  subscripts  again  refer  to  air  and  s~  s   1.      Now,    N„/N 
is  approximately  2,    since  air  is  mostly  diatomic.      The  calculation 
yields  R  =  6.13  Rft.      This   should  be  fairly  accurate  for  helium-3  or 
helium-4.      For  a   1-Mev  proton,    R~  is  2.28  cm,    therefore  R  is    14  0  cm. 
The  range  is  inversely  proportional  to  gas  pressure,    since  density  is 
proportional  to  pressure.      Thus,    at    2.5  atmos  this  proton  range  in 
helium  becomes   5.60  cm,    and  at  5  atmos,    2.80  cm. 

To  find  the  range  of  an  ion  of  other  isotopes  of  hydrogen  (given 

the  relation  for  protons)  we  note  that  of  all  terms  in  the  stopping 

2 
power  formula  only  v     is  different  for  a  given  energy.      Thus,    a 

triton  of  energy  E  experiences  the  same   stopping  power  as  a  proton 

of  energy   l/3  E;  but  it  goes  three  times  as  far   since  it  has  three  times 

the  energy.      The  maximum  triton  energy  with  which  we  will  be  con- 

3 
cerned  in  the  He     disintegration  is  about   1  Mev,    and  we  figure  its 

range  to  be  8.54  cm  in   1  atmos  of  helium.      Since  this  is  less  than  a 

3 

representative  proton  range  from  a  He     disintegration  and  represents 

an  upper  bound  on  the  triton  range,    the  order  of  magnitude  of  the 
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wall-effect  difficulties    can  best  be  estimated  by  considering  the 

[7] 
proton  range.      As    ointed  out  by  Batchelor  in  explaining  his  computer 

program  for  wall-effect  computation,    the  correlation  between  proton 

and  triton  tracks  plus  a  cylindrical  geometry  makes  exact  calculations 

3 
impracticable;  also  the  angular  distribution  of  the  He    (n,  p)T  reaction 

is  not  very  well  known. 

In  integral  bias  counting  involving  the   (n,  p)  reaction,    detection 

efficiency  is  defined  as  the  counting  rate  divided  by  rate  of  emission 

of  protons        When  the  total  range  of  particles  is   small  compared  with 

linear   dimensions  of  a  counter  and  distribution  is  isotropic,    a  value 

for  detection  efficiency  often  given  is 

F(E)  -   1   -  R£/Z  b  -  R£/2  Ls 

where  we  have  an  active  cylindrical  volume  of  radius  b  and  length  L,, 
and  R-p  is  that  portion  of  total  range  of  protons  necessary  to  produce 
a  pulse  above  the  bias  energy  E  of  the  discriminator  .      The  terms 

in  the  expression  for  F(E)  reducing  it  from  unity  represent  the 
magnitude  of  the  wall  effect  and  indicate  the  loss  at  the  lateral  bound- 
ary and  end  face,    respectively. 

The  expression  for  F(E)  can  be  shown  to  be  plausible  as  follows: 
Consider  the  fraction  of  the  volume  of  a  cylinder  of  radius  b  within  a 
distance  R  of  the   side  walls.      This  fraction  is   (2  bR  -  R    )/b    . 
Consider  charged  particles  originating  with  equal  probability  along  a 
radial  line  with  ends  defined  by  the  outer  volume.      If  any  direction  of 
emission  is  equally  probable,    one-half  of  the  particles  will  be  emitted 
with  a  radial  component  toward  the  center  and  not  intercept  the  wall. 
For  that  half  of  the  particles  having  a  positive  radial  component, 
those  originating  on  the  inner  edge  of  the  fractional  volume  have 
practically  zero  probability  of  intercepting  the  walls   (since  only  one 
direction  of  emission  allows  the  particles  to  reach  the  wall),    while 
those  originating  at  the  wall  have   100%  probability  of  reaching  the  wall. 
Assuming  that  this  increase  in  probability  is  linear  as  we  move  out- 
ward toward  the  wall,    we  estimate  that  one -fourth  of  all  the  particles 
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originating  in  the  fractional  volume  reach  the  wall.      Then  the  wall 
effect  is   1/4  time  the  fraction  of  the  volume,    or 

R/2b  -   (R/Zb)2  . 
For  R     <<     2b,    this  is  approximately  R/2b.      Similarly,    the  fractional 
volume  within  a  distance  R  of  the  ends  of  a  cylinder  is  2R/L,    and-- 
again  taking   l/4  of  this  value--we  get  R/2L  as  the  wall  effect  due  to 
the  ends. 

In  our  case  the  condition  of  large  dimensions  is  not  well  ful- 
filled; also  we  assumed  that  the  ionization  from  the  particle  other 

than  a  proton  was  confined  to  a  very  small  region,    which  is  not  so 

3 
true  for  a  He      reaction  as  for  one  involving  heavier  nuclei.      In 

differential  detection  F(E)  is  ideally  the  difference  between  two  step 

functions  with  different  energy  thresholds.      Therefore,    the  deviation 

from  unity  is  a  rough  indication  of  the  magnitude  of  wall  effect  on 

energy  determination  for  a  proton  of  given  energy.      For  a   I-Mev 

proton  in  2.5  atmos  of  helium,    the  geometric  factor  is 

G(E=1  Mev)  =  1  -  R/2b  -  R/2L, 
and  with  2b  =  7.62  cm,      2L  =  38.9  cm,      R  =  5.60  cm,    we  have 
G(E)~0.12  or  the  wall  effect  is  88%.      Doubling  the  pressure  would 
reduce  the  wall  effect  to  44%,    and  at   10  atmospheres  to  22%. 

For  counting  neutrons  in  ordinary  helium  we  are  concerned 
with  the  range  of  the  alpha  recoil.      The  maximum  energy  of  the  alpha 
particle  is   16/25  E    ,    as  shown  in  Chapter  VIII.      If  we  consider;    the 
average  energy  of  the  Pu-Be  neutron  spectrum  to  be  4  Mev,    this 
corresponds  to  a  maximum  alpha  energy  of  2.56  Mev,    or  an  average 
E     of  about  half  this  amount.      Figuring  the  range  of  this  alpha 
particle  in  2.5  atmospheres  of  helium  by  our  previous  method,    we 
get   1.38  cm.      The  wall  effect  is  22%  for  this  range.      Experimentally, 
a  coincidence  counting  rate  of  about  16%  of  the  total  counting  rate  in 
the  main  counter  was  measured.      This  is  a  very  rough  check,    but  it 
at  least  appears  that  the  ring  counter  is  capable  of  reducing  the  wall 
effect. 
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ELECTRON  RANGE 

To  obtain  the  range  and  energy  loss  of  electrons  in  estimating  the 

effect  of  tritium  decays  or  background  gamma  counts     one  can  take  the 

[41  1 
t  tabulated  value  of  stopping  power  for  helium  as  a  function 

2 
of  electron  energy,    given  in  units  of  Mev  per  g/cm    .      For  example     at 


nearest  tabulated  value  ot  stopping  p 

of  electron  ener 
10  kev.    we  have 


-  dE/dx  =  22.5  Mev/g/cm 

If  we  take    A  x  as  the  average  path  in  the  counter,     1,59  cm.    times  the 

-4      /        2 
density  of  helium  at  2.5  atmospheres,    Ax  =  7.53  •     10        g/cm    .      Then 

AE  is  greater  than   16  9  kev  (the  product  of  -  dE/dx         Ax),    because 
the   stopping  power  increases  with  decreasing  energy.      For  higher- 
energy  electrons  it  is  fairly  accurate  to  use  Feather's  Rule  or 
Flammersfeld's  formula  *•  to  find  x  and  then  multiply  Ax/x  times 

the  energy  of  the  electron. 
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CHAPTER  VII 

NEUTRON  GENERATOR 

-Eor  testing  purposes  one  desires  a  monoenergetic  neutron 
source.      The  complex  spectrum  of  (a,  n)  sources  has  been  mentioned; 
the  other  class  of  neutron  sources  using  radioactive  nuclei  are  the 
photoneutron  or   (y,  n)  sources.      These  use  deuterium  or  beryllium 
because  other  nuclei  have  thresholds  above  6  Mev.      To  get  high  enough 
intensity  the  gamma  source  is  surrounded  so  that  there  is  an  intrinsic 
energy  spread  due  to  different  angles  between  neutron  emission  and 

incident  gamma  ray.      This  is  small;  for  example,    it  is  3.37.%    for  a 

24  9  [201 

Na        +  Be      source1         .      However,    this  is  not  the  main  cause  of 

energy  spread.      Owing  to  the  large  quantity  of  beryllium  there  is 

considerable  neutron  scattering  in  the  source,    and  gamma-ray  loss 

of  energy  by  Compton  scattering. 

It  has  been  found  possible  to  construct  neutron  generators 

3 
utilizing  the  deuter on-deuteron  reaction,    D(d:,  n)He    ,    with  useful 

neutron  intensities  without  elaborate  high-voltage  equipment;  for 

[42l 
example,    the  generator  of  Zinn  and  Seely  was  built  in  1937. 

Voltages  on  the  order  of  100   kv  can  give  useful  yields  even  though  the 
reaction  is  not  very  efficient  at  this  energy.      With  this  low  voltage  we 
also  have  the  advantage  of  almost  monoenergetic  and  gamma-free 
neutrons.      For  a  thick  target  the  deuteron  energy  at  the  time  of  the 
reaction  can  vary  from  a  maximum  to  zero.      In  the  generator  used 
in  this  experiment  the  maximum  deuteron  energy  is   120  kev.      Thus, 
using  equations  developed  in  Chapter  VIII  and  Q  =  3.28  Mev,    we  find 
that  in  the  forward  direction  the  neutron  energy  can  vary  from  2.46 
to  2.91  Mev.      At  a  laboratory- syfetem  angle  of  90  degrees  the  energy 
spread  is  from  2.46  to  2.49  Mev,    or  little  more  than   1%.      The  energies 
are  weighted  toward  higher  values  because  yield  increases  rapidly 
with  deuteron  energy.      This  neutron-energy  range  is  convenient  for 
checking  the  operation  of  the  counter  because  the  maximum  recoil- 
alpha  energy  is  about  the  same  as  the  maximum   (E      f  Q)  energy  we 
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wish  to  measure  in  helium- 3,    and  about  half  the  recoils  are  above  the 
lowest  sensitivity  needed. 

PULSED- NEUTRON  GENERATOR 

The  neutron  generator  that  is  being  put  into  operation  provides  a 

pulsed  source   of  neutrons  and  is  similar  to  the  one  described  by 

[43] 

Ruby  except  for  the  method  of  mounting     the  accelerater   system  and 

the  vacuum  equipment.      The  accelerator   system  is  described  in 
Reference  43;  briefly,    the  deuterons  are  generated  in  a  pulsed  Phillips 
Ion  Gauge  type  of  ion  source   supplied  with  deuterium  gas  through  a 
palladium  leak  and  accelerated  to  impinge  upon  a  circular  target  of 
l-l/8-inch  diameter.      The  targets  are  prepared  by  occluding  deuterium 
in  titanium.    This  type  of  target  allows  simpler  equipment  than  one 
with  heavy  ice,    by  requiring  no  special  cooling.      It  is  planned  to  try  a 

deuterated  paraffin  target  to  see  if  the  yield  is  increased.      The 

1Z  ,  13 

objection  of  non-monoenergetic  neutrons  due  to  the  C       (d,  n)N        re- 
action is  not  pertinent  here  because  we  are  below  the  reaction 
threshold.      The  generator  was  designed  so  that  adjustment  of  param- 
eters would  not  be  critical;  the  yield  is  largely  determined  by  the 
accelerating  voltage  and  the  condition  of  the  target.      (The  reaction 
cross  section  decreases  with  deuteron  energy,    and  the  targets  have  a 
tendency  to  form  a  carbonous  surface  after  a  few  hours'  use.  ) 

The  pulsed  operation  is  particularly  convenient  for  experiments 
that  require  a  reference  time,    such  as  in  neutron-diffusion  experiments; 
single  pulsing  has  been  used  for  cloud  chamber  experiments.      In  normal 
operation  the  ion  source  is  pulsed  at  60  cycles  per  second  with  200- 

usec  pulses.      In  this  manner  we  can  get  an  integrated  yield,    Q    ,    of 

6  S 

more  than  10     neutrons  per   sec.      This  is  a  duty  cycle  of  0.01Z,    or 

during  a  pulse  the  instantaneous  neutron  yield  is  Q   /0.012  neutrons 

per  second.      For  a  cloud  chamber  experiment  we  used  an  external 

oscillator  option  of  400  cycles  per   second  and  increased  the  width  of 

the  pulse  by  an  adjustment  in  the  arc  pulser  to  500  u-sec.      This  is  a 

duty  cycle  of  0.2.      This  mode  of  operation  was  gated  on,    when  it  was 
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desired  to  take  a  picture,    for   100  milliseconds,    giving  a  total  pulse-on 
time   of  20  msec  per  gate.      Then  the  neutron  yield  per  gate  should  be 
Q   /0.012  times  20  msec,    or  more  than  10     neutrons  in  each  burst, 
assuming  the  same  effectiveness  of  the  ion  source  when  pulsed  at  400 
cycles  per   sec.    as  at  60  cycles.      It  would  be  more  realistic  to  calibrate 
at  400-cycles   steady  operation,    but  this  increases  the  risk  of  damage 
to  the  ion  source. 


NEUTRON  YIELD 

In  order  to  measure  the  integrated  yield  one  should  determine 
if  the  amsotropy  of  the  yield  is  significant.      It  has  been  found  that  in 

the  energy  region  below  0.5  Mev,    the  neutron  yield  as  a  function  of 

[44] 
angle  in  center- of-mass  coordinates  can  be  represented1  by 

N(0)  =  A  (1  +  B  cos2  6), 
where  A  is  a  constant  and  B  is  a  number  that  increases  with  deuteron 

energy.      Thus,    the  yield  is  smallest  at  6  ~  90  degrees.      The  quantity 

[45l 
B  has  been  determined  experimentally  down  to  500  kev,    and  the 

curve  extrapolates  fairly  well  into  a  theoretical  determination  at 

3 
low    energies  based  on  the  D(d,  p)  H     reaction.      It  is  believed  the 

angular  asymmetry  is  the  same  for  the  companion  reactions  at 

corresponding  voltages.      We  estimate  B  as  0.6  at   120  kev.      This 

means  the  neutron  flux  at  6  =  0  should  be  (1  +  B)  times  the  neutron 

flux  at  9  -  90  degrees.      Or,    performing  an  integration,    we  find  that 

if  Q     is  determined  by  measurements  at  6  -  0,    we  should  divide  this 

by  (1   -KB/2+B))  for  the  true  yield  or   1.23  for  B  =  0.6.      This  is  a  high 

estimate  of  the  anisotropy,    because  not  all  the  deuterons  have  the 

maximum  energy.      With  these  small  deuteron  energies,    we  expect 

that  the  laboratory- system  angles  will  not  be  much  different  from  the 

center-of-mass  angles.      We  need  not  outline  the  calculations  here 

(a  description  of  the  two  coordinate      systems  and  some  convenient 

equations  may  be  found  in  Reference  47),    but  for  6  =  90  degrees, 

we  find  the  lab  angle  6       =  85.6  degrees  for   120  kev,    and  for    '0 *  -90 

degrees,    6  =  94.4  degrees. 
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The  generator  yield  has  been  monitored  with  a  BF.  proportional 
counter   surrounded  by  2.5  inches  of  lucite  moderator  and  a  thin  sheet 
of  cadmium.      Care  must  be  taken  in  its  calibration  since  this 
arrangement  does  not  give  equal  sensitivity  over  a  wide  range  of 
neutron  energies.      For  a  particular  counter  we  found  a  sensitivity  of 
0.29  counts  per  neutron  per  cm     when  a  Pu-Be  source  was  used,    and 
0.37  counts  per  neutron  per  cm     when  a  mock  fission  source  with  a 
smaller  energy  range  of  neutrons  was  used.      However,    it  can  be 
calibrated  on  the   D-d  neutrons  by  comparison  with  a  calibrated 
Hanson  and  McKibben  type  long  counter.      The  lucite  embedded 

counter  is  easier  to  place  and  adds  less   scattering  objects  to  the 
experimental  arrangement. 
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CHAPTER  VIII 


NUCLEAR  REACTION  DYNAMICS 


©— >V- 


Consider  the  general  case    sketched    above  of  a  particle  of  mass 
M,   and  velocity  v.    colliding  with  a  target  nucleus  of  mass  M?  assumed 
at  rest,      After  the  reaction  the  resulting  particle  of  mass  M_  flies  off 
with  velocity  v_  at  angle  0T  .       The  nucleus  of  mass  M.  has  velocity 
v  .  at  angle  <j>T  .      We  will  use  nonrelativistic  dynamics,    mass  numbers 
for  masses,    and  laboratory  coordinates. 

First,    write  two  equations  for  conservation  of  momentum  for 
components  parallel  and  perpendicular  to  v.: 

M   v     =  M,v       cos    4>L  +  M3V3  cos  0     ,  (VIII-  1) 

M3v3  sin  0L  =  M4v4  sin  <|>       .  (VIII-  2) 

By  transposing  the  last  term  of    Eq.    (1)  to  the  other   side, 
squaring  both  equations  and  adding,    we  can  eliminate  <}>  T  .      Then 
replacing  velocities  by  use  of  E  =    -|  M  v        for  energy  and  rearranging, 


we  obtain 


M,  M. 


E4  =     tC      El   +    S"   E3  "    NT    (E1E3M1M3)i  cos    6^ 
4  4  4 
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We  can  write  another  equation  for  E     from  the  definition  of  the 
Q  of  the  reaction,      Having  assumed  E?  -  0,    we  write 

E4  =  0  +  Ej   -E3. 

3 
For  the  He      (n,  p)T  reaction,    using  Q  =  0.77  Mev  and  similar 

symbolism,    we  write  two  equations  for  the  triton  energy: 


ET  =  0.77  +  E 
1  n 


With  these  two  equations  we  can  solve  for  the  proton  and 
triton  energies  for  given  values  of  E     and  61  . 

For  En  =  0,    the  proton  receives  75%  of  the  Q  energy  and  the 
triton  receives  the  remaining  25%.      That  is, 

E     =  0.577  Mev, 

P 
ET  =  0.193  Mev  . 


For  E     -   1   Mev,    a  tabulation  is  given  below. 


(degrees ) 


Proton 


E  %  of 

P 
(Mev)  (En  +  Q) 


Triton 


T 

(Mev) 


%  of 
(En  +  Q) 


0 

1.74 

98.3 

0,03 

1.7 

45 

1.  52 

85.8 

0.25 

14.2 

90 

1.08 

61.0 

0.69 

39.0 

135 

0.77 

43.5 

1.00 

56.5 

180 

0.67 

37.8 

1.10 

62.2 

It  is  often  desirable  to  have  a  convenient  expression  for  the 
energy  of  the  ejected  particle  at  6  ,    =  0  and  at  90  degrees  as  a  function 
of  the  bombarding  particle  energy.      For  the  forward  direction  the  two 
equations  become 
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M 


E.   = 


M. 


4       MJE1+    M7E3-i5-     (E1E3M1MJ)>     , 

4  4  4 


E4  =  Q  +  E1-E3 


A  solution  for  E     yields 


E3  =  A2/2A3  +  \   [(A2/A3)2  -  (4A^/A3)]  »  , 


where 


Aj  =  M4Q  +  Ej  (M4  -  Mx)  , 

A2  =  2AX  (M3  +  M4)  +  4  Ej  M1  M3  , 


A3  =  (M3  +  M4)' 


For  the  He      (n,  p)  T  reaction  this  gives 
EP=    I    Q    +   fEn<1  +   it9+     TT^-]"> 


For  the   D(d,  n)  He     reaction, 

E     =-Q+-E     (l+-[3+    6Q     12   ) 
n       4  U  +    2        d  u  +    2  L^       E,      J      '  ' 


For     6,    -  90  degrees,    the  two  equations  are 

M  M 

E4  =    M7  El  +   MT     E3  ' 
4  4 


E4=  Q  +  Ej  -   E3 


Solving,     we  obtain 


E3  = 


M4/(M3  +  M4) 


Q  +  [(M4-M1)/  (M4  +  M3)]     E1 
For  the  two  reactions  above,    at  90  degrees, 
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3  1 

E     =  i.  Q  +4    E 
p       4  2        n 


3  1 

E      =  ±Q  +i-    E, 
n       4  4       d 


For  the  case  of  elastic  scattering  we   set  Q  =  0  and  M_   =  M 
therefore  the  two  equations  become 

MM  i 

4       M.         1         M.        3  M.  13        1  L 

4  4  4 

E4  =  El  "  E3  ■ 

We  are  interested  in  the  energy  given  to  the  recoil  nucleus, 
therefore  we  eliminate  E-  from  the  first  equation  by  using  the  second. 
Then  we  set  0T    =   180  degrees  for  maximum  E  .  corresponding  to  a 
backward  recoil.     Solving  for  E  .(max),    we  have 

4M,M4 

E      (max)  =    ?      E       ; 

(M1+M4r 

For  neutron  scattering  in  helium- 3, 
ERe3  (max)  =    §    Eft  ; 

and  for  neutron  scattering  in  helium-4, 

E      (max)  =    -=^    E 
a  25        n 
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